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Foreword 


When the ERTS px-’oject was officially announced by NASA it was 
received with great interewSt in Peru since it represented a way to 
develops Remote Sensing Techniques, that could substantially acoeler 
ate and improve the study of our natural resources, without the 
higher initial investment necessary to acquire airborne multiespec-- 
Iral sensors. 


The Geophysical Institute of Peru (IGP) and the National Office 
for Evaluation of Natural Resources INERN) joined their efforts for 
Remote Sensing development and with the help of the National Agra-- 
rian University prepared a proposal to pax'ticipate in the ERTS pro- 
ject and presented it to NASA in 1971. 

Since then an Intense program of activities was developed and 
in spite of initial difficulties it has evolved toward a succes- 
fui program. 

Previous to the launching of the first ERTS satellite, techni- 
cal groups were organized to work on the ERTS imagery that x^ould be 
arriving, all the basic information that was available on the Santa 
River Basin was collected and at the IGP, Mr, Eduardo Seminario un- 
der the direction of Dr. Medardo Molina from UNA made a Hydrologic 
Study of the Santa River Basin to support the studies with Landaat 
Imagery. 

In order to promote the national interest on remote sensing and 
increase support to this project, talks were given at Universities , 


national meetings and government organizations. These efforts 
culminated with a Seminar on Remote Sensing given at UNA in 
Ajjril 1972. 

in July 1972, an agreement between ONERN and IGF was 
signed to work the development of Remote Sensing in Peru and to 
-eek national or foreign finantial aid to support it. 

ER’ S-1 wa.-: suooesfully launched in August 1972, after 
a number of postpcnements , and finally the first images started 
to arrive during the first half of 1973, 

The images v^ei’e taken over our observation site on the 
Amazon Ji.ingle. At that time oil exploration activities in the 
jungle were booming and a prioritary interest was given by the 
government to that region of the country. It was then decided 
that the efforts of the teams at IGF, ONERN and UNA will concen 
trate their work on studies about the jungle areas covered by 
ERTS-1. 


Mr. Eduardo Ceminario and Dr. Molina made a study of 
Meanders over rivers in the Amazon Jungle. By the end of 1973 
Mr. Seminario v^as leaving to Europe to benefit from a scholar-- 
ship and Mr. Molina took a leave of absence in early 1974 to 
woi'k for UNESCO. 

At IGF an infraestructure was implemented to deal with 
the imagery coming frem NASA, to enlarge the black and white 70 
mm. and to produce 3S mm color composites. 
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The imagery received from NASA was given broad diffusion 
being the most significant applications the use of 250,000 en- 
largements as a help for a photomosaic work at Servicio Nacional 
de Aerofotografla (SAN), the use of 1:250,000 enlargements by 
Belco Petroleum Co. in their planning of an oil pipeline in the 
NortheiTi peruvian jungle. The use of 1:500,000 enlargements by 
Minero Peru in their research work. 

In 1974 under the auspices of the National Institute of 
planning a proposal by ONERN and IGP was presented to the Cana- 
dian International Development Agency (CIDA) to finance a pilot 
Remote Sensing project. The preliminary part of this project 

is already underway and the main part is expected to start 

this March with an study based on Landsat 1 and 2 imagery. 

This report presents the work that was performed at 
IGP, ONERN and UNA during the period of 1972 tlirot^h 1975 with - 
the Landsat 1 (ERTS-1) imagery. 

During the course of this work it was decided that ONERN 
will concentrate their efforts in the development of photointer — 
pretation techniques for multiespectral imagery as a natural evo- 
lution of their work on the evaluation of natural resources of 
the country. IGP will work mainly in computerized methods of di^. 
ta handling an analysis. UNA concentrated their efforts in hy- 
drologic studies. 
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In the last eight yearti a great deal of effort have been 
spent in oil exploration activi.tief;; in the peruvian iMnazonic jun- 
gle. At this time a pipeline that will bring, the oil from the 
northern fields of the jungle to the coast is under rrjn struct ion. 
With the of the oil crisis still present, is not diffi- 
cuic to understand the interest of Peru in this vast region of 
the couMi’y. 

With the urgency to develop our capability to produce 
oil came ttie necessity to obtain basic information about the jun- 
gle. The first images from ERTS-1 on that area showed that a 
great deal of ijnportant information could be obtained from them. 
The Geophysical Institute of Pe.u (IGP) , the National Office for 
Evaluation of Natural Resources (ONERN) and the Agraria Univei'si. 
ty (.UNA) stax’ted working with those ERTS-1 images', ONERN and 
IGP v/orked to develop olassifical ion methods for* multlespectral - 
imagery and UNA in the study of meanders in the Amazon Basin. 

The first task was to select imagery with good quality, 
low cloud cover and with a largo variety of intei’Cst i.ng featui'es. 
The image selected was the one taken over the region close to the 
cities of Contamana and Tarapoto. 

To develop classification schemes, photographs from 
airplanes were used along vtfith enlargements to 1:2S0,000 scale 
and color composites of 1:1,000,000 scale. In the course of 
this v^^ork It was found that significant detail was being lost in 
the photo interpretation process. The problem was mainly due to 
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i:ric; of the following causes: the photo interpreter found difficu]^ 

!;iea to sitfitoh from one channel to other to identify a ‘ feature, 
the details were too small on the 1:1,000,000 color composites or 
too faint in the enlargements, the photo interpreter got tired and 
started making errors. 

It was then decided that computerized technj.ques should 
studied as a way to solve some of these problGms and to aceel^ 
ral v: the interpretation process. The initial results at ICP were 
pronii ing and then a moi'e systematic effort was started in this 
field. It was decided that ONEiyvl will continue working with pho- 
to interpretation methods and IGF will develop computer interpre- 
l at ton techniques, bot:h works were carried out in paral.u i. 

The purpose of this report is to present what was accom- 
plished at the IGP to develop computer techniques to process eoni- 
(iuter compatible tapes (CCT) from Landsat Satellites, 

This work concentrated in the multiespectral analysis of 
CCTa, the problem of geometric correction and application of tech 
niques to define geometric features and temporal variations was 
left for future development. 

V\fhen this work .started a Uatacraft 602h comparer was a- 
vallable at IGP, the memory capacity was 32 K bytes and the tele- 
printer was too alow for practical use which made ihe progress too 
difficult. 

Then, it was decided to use the IBM 3GO/40 computer at 


ORIGINAL PAGE IS 
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the National Unlverrjity of Bngineering, but it was found that a 
large computer center is not particularly suitable for the type oi' 
work that we were interested in. Finally a Hewlett Packard 2100 
computer systeirs installed by Stanford Research Institute at the 
Ancon Observatory of IGP became available and then it was possi-- 
ble to satisfy most of the requirements for our work and the pro- 
grams for machine processing of CCTs reached a more satisfying 
statu:.. 


Report gummary 

In chapters HI and IV the general problem of machine 
classification is introduced followed by a brief description of 
CCTs in chapter V. 

The computer techniques that were developed dui'ing this 
work are presented in chapter VI and their application to the 
study of jungle areas are described in chapter VII. 

The conclusions that were reached are presented on chap- 
ter VIII, and finally chapter IX describes our plans for future -■ 
work. 

III. The Problem of Machine Classification of Remotely Sensed 

Data 

Fig. 1 ilustrates the general problem of machine classi- 
fication. The source illuminates the terrain, the sensor recelv'cs 
the reflected radiation and at a particular time generates a so i 
of numbers, patterns, which are a measurement of the received i\i- 




SOURCE PATTERN 



TERRAIN 

(q) 



FIGURE I 

CLASSIFICATION OF REMOTELY SENSED DATA 

INTO N CLASSES 


I 

{ 





diation. 
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The feature selector extracts those numbers that are sig 
niflqant to define or characterize the surface under observation. 
This ncxv set can be represented by a vector generally known as the 
feature vector. 

The purpose of the classifier is to assign each input pat 
tern to one of the possible pattern classes. The classifier malces 
this assignment on the basis of the feature vector obtained from 
the input patterns and also based on past experience. 

The development of a classifier has been the main objec — 
tive of this work and the feature extractor has been only consi- 
dered as a particular c.i.se of the classifier. 

Xno problem of classification is basically one of parti- 
tioning of the space defined by the feature vectors in regions, 
one region for each class. The partition is performed in such a 
v^;ay that the probability of missclasif ication is minimized. Ex- 
tense literature exists on the various methods to solve this pr^ 

At present there are a number of methods for machine cla^ 
4 6 

sification ' , they are based on the use of the covariance ‘matrix, 
in the measurement of probability distributions (parametric and 
non parametric) clustering techniques and other. The method cho- 
sen will depend in general on the problem itself and on the avai- 


*Super?.cripts indicate references at the -end of report. 
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3 able digital equipment. 

The main purpose of this work was to develop computer tech 
niques to analyse Landsat-1 imagery over the Peruvian Amazonic Jun- 
gle close to Ucayali River, one of the Amazon tributaries. The 
structure oi the multiespectral data obtained was very complex, ob 
iously non gaussian, ground thruth is almost nonexistent and dif- 
f i.'ult to obtain because of the inaccesibility of the region. The 
digital equipment used in the problem was based on a Hewlett Pa- 
ckard 2100 computer \\;ith 32 k bytes of memory, two tape units and 
a 600 character/min teleprinter. These facts led us to choose 
clustering techniques as the base for our classification algorithms. 

The Cl u stering Method 

Clustering techniques are based on the determination of 
Natural groupings in a set of data. The space defined by the fea- 
ture vectors is structured and there are regions where the density 
of points is .greater than the surrounding space, clustering me- 
thods permit to define volumes that will enclose those regions. 

Intuitively, one can assume that data points, or feature 
vectors, corresponding to objects on the terrain that are similar 
sliould come close in the feature space and hopefully forming clu^ 
tens. 

It is fundamental to the clustering metnod the definition 
of a metric to determine closeness, or similarity, in the feature 
space. From the variety of similarity measures'^, the Euclidean 
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distance between feature vectors was used in this work and close — 
ness was define in base of the nearest neighbor^^ criterion. 

V, The Computer Compatible Tape fCCTl 

As a new tool for the study of remotely sensed land sur- 
faces come the CCTs^ generated from data collected by the multi- 
espectral scanners (MSS) on board of Landsat Satellites. 

The MSS covers the visible and near infrared spectrum - 
with four channels: 

Channel 

4 0,5 - 0.6 microns 

5 0.6 - 0.7 '* 

e 0.7 - 0.8 ” 

7 0.8 - 1.1 

The MSS scans cross track swaths 185 km wide at a nomi- 
nal altitude of 494 nautical miles. It consist of six detectors 
for each of the four channels, video outputs from each detector 
are sampled, commutated and multiplexed. The commutated samples 
ai’C encoded and transmitted to ground-base receiving sites. The 
receiving stations compile the raw data on video tapes and trans. 
mit these tapes to the NASA Data Pracessing Facility (NDPF) at 
the Goddard Space Flight Center, Greenbelt Maryland. The NDPF 
corrects, calibrates and formats the raw MSS data and converts 
it to binary form on Computer Compatible Tapes (CCT) . 


A 185 km square ground scene on the CCT is a final pro- 
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tlucL of the MS.S. fig, 2 illustrates a typical scene correspond- 
ing to a pass ul r»and.sat~l Satellite over the Ucayali and Hua- 
llaga River in the Central Amazon Jungle of Peru, 

A scene is made up of parallel scan lines each contain- 
ing a large number of video xioints or picture elements (Pixels). 

The NJ)Pl' recurds a complete ground scene on four sepa- 
rate CCTs, Each CUT contains image data for one 46. 25-185 Kni 
strip. The CCT contains data that has been radiometrically ca- 
librated^ but not geometrically corrected for skew due to earth 
rotation. CCTs come in 800 bpi and 1600 bpi densities of tapes. 

In this work 1600 bpi CCT corresponding to the image 
:>hov,’n in figure 2 v;as used. 


VI. The bcveloomen t of a Computer Techniques to Analyse Re- 

motely Sensed Data Recorded on CCTs 

As a framework for this work the following basic requi- 
rements were defined: 

(!) Clustering techniques will be used, 

0 Tlie progi*ains will be interactive so that they permit, 
as much a.s possible within the limitations of our 
equipment, the intervention of the specialist during 
the analysis. 

0 The programs should use the four channels from the 
CCTs but will have, as an option, the use of less num 



FIO 2 TYPICAL MSS SCENE FROM LANDSAT f 


CHANNEL 5 
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her of channels. 

vO The presentation of results will be, either on te- 
leprinter or CRT, 

>1 Although they will be developed with a particular 
pr'oolem in mind, they should be sufficiently gene- 
ral that they can be used in other applications 
tor . 

■|iu; rogra'.is should use, as much as possible , stan^ 
;u'd languages, like ALGOL or FORTRAN IV. 

Fig. 3 illustrates the equipment used in this work. 

The group of programs developed to process CCTs was named 
INT1«RPRL‘1‘ ami consist of the following; 

hiubroutines for data manipulations from CCTs, 
written in ASSEMBLER language callable from ALGOL 
and FORTRAN IV programs. 

Programs ALMAP I and ALMAP II display on the tele- 
printer and oscilloscope respectively the video .in 
formation, levels of reflectance, from. one of the 
channels stored on CCTs. The display is in the form 
of a selected area from a scene. The language used 
is ALGOL, 

These programs have an option that permits to li-- 
mit the range of the video information so that, 

only those areas on the terrain that produce rc-- 
flcxions JeveJs within given limits will be dis-- 
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played. Also the desired range, between the limits 
can be divided into slices of equal width and, 
displayed with different symbols on the teleprinter 
or different gray levels on the Oscilloscope, this 
permits to identify features on the terrain with 
levels on the video information. 

It v;as found desirable to be able to vary the scale 
on the display, and therefore ALMAP I and II were 
provided with an option that permits to produce di^ 
plays at a prescribed scale. 

The programs are interactive through a continous 
interchange of messages between the machine and 
the operator. Fig. 4, shows a sample of ALMAF I 
messages. Fig. 5 and 6 show typical displays from 
ALMAP I and II respectively. 

Programs HISTA and HISTB display, on the oscillos- 
cope and teleprinter respectively, two dimensional 
histograms of any pair of MSS channels. The value 
of histogram elements are indicated by means of 
gray levels on the oscilloscope or symbols on the 
teleprinter. These programs are interactive and 
have options that permit the user to display a pre^ 
cribed range of levels or to slice out the more 
intense peal<s so that to enhance the fine structure 
of the histograms. 


ALGOL MODIFIED MAP PROGRAM (AWCON SEPTEMBER 75) 

UNIT, TAPE?, RECORD LCNGHT? , RECORD POSITION? 

9.3296,-1 

TITLE? 

NORTE OE CONTAMANA 
CHANNEL? 

6 

INITIAL LINE?, FINAL LINE? , INITIAL POINT?, PINAL POINT? 
25fl,320.2Ofl,260 

LINE INTERVAL?, POINT INTERVAL? 

1,1 

MINIMUM LEVEL? , MAXIMUM LEVEL? 

0,33 

LEVEL INTERVAL? 

3 

CHARACTERS? 

SI 

000300300111 
123456789012 
abcdefghi jkl 
NORTE DB CONTAMANA 
CHANNEL: 6 

PROM TO CHAR 
0 2a 
3 5 t» 

6 8c 
9 11 d 
12 14 e 
1*^ 17 f 
I.'- 2.1 n 

21 23 h 
24 26 i 
27 29 j 
30 32 k 
33 33 1 


FIG. 4 - Program A1J4AP I. Example of messages 
interchanged between the user and the 
computer. ‘ 
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f'ig. 5 Map of Ucayali riv6r. Program ALMAP I, characters 
to different gray levels in channel 6. Vertical ' 
indicate scan I Ines- 
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FIG. 6 - MAP OF UCAYALI RIVER WEST OF CONTAMANA 


PROGRAM " ALMAP II" 
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® Program CLUST is the main one, it permits to de-- 
fine clust'irs in the feature space and has been de. 
veloped based on the technique named ISODATA (Iter, 
ative Self Organizing Data Analisys) , Reference 5. 

ISODATA is a collection of iterative techniques - 
that permit to find the principal modes, or sub- 
clusters, in the data by grouping those patterns 
that are close together, and representing them by 
a hypersphere centered on the average of the pat- 
terns. 

By means of these techniques a' cluster is repre — 
sented by a number of subclusters in such a way 
that the variance of each subcluster is kept- small 
er than a parameter defined by the user. In each 
iteration the sLibclusters are' split or lump toge- 
'ther depending on a set of control parameters. Fig 
7 and Table I briefly describe ISODATA. 

CLUST has been implemented in- ALGOL with a modular 
structure, it consists of a set of subroutines 
that perform each of the functions described in 
Fig. 7. CLUST is interactive and permits the oper 
ator to redefine all the control parameters after 
a given number of iterations. 

The results obtained by CLUST are a set of subcluj^ 
ters which are printed in form of -a table specify- 









ABLE 


□DATA Main I unctions as Specif led in Fict. 7 


Step 


Function Description 


Initial .subcluster centers are fed to the 
t.'omputer anu the control process parameters 
are specified. 


Input patterns are allocated to the subclus- 
ters using a minimum distance criterion. 

Ti^e following variables are computed: 

Average vector for each subcluster, 
Standard deviation for each sub — 


olur;l-"r, 


- Average distance of the patterns to -their 


subcluster center 


Numbf:r of patterns for each subcluster 


Discard tho;;c suliclusters whose number of pat- 
terns is Icsq than the parameter On. 


Take this way if; 

Last Iteration 
Or even iteration 

Oi tumber of subclusters greater than 
two times the number of subclusters (2 N) 
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T A B L E I 


(Continuation) 


Step 


Function Description 


specified at the beginning of the process. 

Talce this way if : 

-- Odd iteration 

- Or number subclusters less than half the num- 
ber of subolusters (0. 5 N ) specified at the 
beginning of the process. 


Lump together two subclusters if the distance between 
their average vectoi's is less than 0c. 


Split a subcluster in two if in any dimension the 
standard deviation is greater than 0e, 
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ing their main characteristics such as feature 
space coordinates, population, standard devia- 
tion and radius, fig- 8 shows a typical prin^ 
out from CLUST. 

. Program CLASIF written in ALGOL performs the 
classification based on the subclusters de- 
fined by CLUST and selected through an opti- 
mization process that will be described later. 

A pattei'n is assigned to a subcluster if its 
Euclidean distance to that subcluster is the 
minimum distance to any subcluster and it is 
smaller than a maximir -adius defined during 
the optimization process. The result of this 
classification can be displayed by the tele- 
printer or the oscilloscope. 

A modification was made to the teleprinter 
symbols so that to produce a gray level presen 
tation on the teleprinter. 

The program is interactive and permits to 
change the area to classify and the subclusters 
parameters»Fig.'9 shows a sample of CLASIF mes- 
sages. 

With these set of programs it was possible to develope a me- 
thodology to perform the computer analysis of CCTs. Figures 

10, 11 and 12 describe the process of analysis in which the user 
plays a central role. 
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Fig. 9 Typical messages betvveen the user and 
the computer. Program CLAS^F, 




TABLE II 

Symbols used on fXo\v charts 


HP 2100 


HP 2100 Computer 
32 K Bytes 



Magnetic Tape unit 
Loaded with a CCT 



1200 Diablo Printer 



1200 Diablo Keyboard 



J 


Tektronix Oscilloscope 
Equipped with a Polaroid . 
Camera 



Remex Fast Paper Tape Reader 
and Punch 


C.specialist usin>f 
INTERPRET 


HEPHv. . 
PRIGINA. 











STATISTIC FROM 
PREVIOUS STEP 


a 


COORDINATES 
OF SAMPLES 













CLUSTERS 

POINTS 


HP 2100 
PROGRAM 
CLASS! F!C 


THEMATIC 
. MAP 



FIG 12 FLOW CHART FOR COMPUTER CLASSIFICATION 
CONTINUATION FROM PREVIOUS PICTURE 
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Oriee the- problem has been de;iined the process of anaLysis 
•starts with the selection of. samples.. The us-er helped by all the 
information available, aerophotography and ..ground thruth, selects 
areas , with hnovm characteristics by displaying the. Landsat .infor- 
mation ohvthe teleprinter, by. means of ALMAP I. 

The user determines the. coordinates-' of the samples and ask 
the eomputer to. punch tliem on paper -tape. , With the help of HISTA 
or HISTB the user display on the teleprinter, or on a Polaroid - 
picture, two dimentional- histograms obtained -for the sample areas, 
later thase-,hxstograms vjlll be used to . check the results of - 
GLUST. 

The second step of the process is to determine subclusters 
from the sample data by means of CLUST. Once the subclusters are 
determined they are mapped on two dimentional histograms and com- 
pared with the histograms determined in the previous step. 

, The sample patterns that are fed to CLUST are not comple- 
tely "pure", that is, they contain data that do not belong to a 
known classification and therefore those patterns will be rejected 
during tliis step. 

Next the sub clusters are passed through an optimization pro. 
cess by means of Che program GLASI'F, Each subeluster is tested 
using the classification algorithm of CLASIF and the resulting maps 
are displayed on the bseilloscope , and on the teleprinter. Data 

from a. loioim area, , different from that of the sample data, is used 

■ . . . ; 

•for this process. , The radius of each suhnluster is reduced to mi- 
nimize the number of missclassifications or increase to maximize 


: 2-2B', 

correct olasstficatioris. 

The last step is the analysis of the area of interest i the 
classification is performed using GLASIF and the resulting maps 
■ are eitlier displayed on the teleprinter or on oscilloscope as il- 
lustrated on figure 12. If the maps are not satisfactory then new 
sampling areas are defined and the process is repeated again but 
only limited to those classes that were not being correctly cla^ 
sified. The process usually, does not require:. more. .than two itera- 
tions. 

VII. Applications of IMTERERET to the Study of Amazonian Jungle 

The jungle of Peru is scarcely developed and its present 
conditions, in most of the amazon basin, is the same as it was 
millions of years agp^^; ;:.; I^^ evolved to adapt it 

self to the environment (soil, geomorphology, weather, etc.) and 
has reached the condition of climax forest in most of the Jungle. 
The rivers are the m^in dinamic agents and the wavelike evolution 
of the meanders create typical patterns with oxbow lakes and for- 
mer river beds covered with heryapeous vegetation. 

After discussions with professor Jorge Malleaux from the - 
Forest Department of the Agrarian National University the area 
close' to the town of Contamana on the Ucayali river was chosen , 
zone I of Fig. 13. This region was attractive for this work be- 
cause of the variety of. forest types that presented. 

At first, v^hile INTERPRET was still in its preliminary form 
it was used for an identification of drainage patterns, lakes, man 
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FIG. 13 Scene from Landsat 1 chosen for this study zone I 
and zone II are the area covered. 
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made clearirigs of tlie jungle or places where herbaceous vegtjtation 
has started to grpw, swamps and wet areas. 

In October 197 under the auspices of the National Aero~ 
pihotographic Seryic (SAN) a trijp to zone I of Fig. 13 was made 
by persbnhei of ON arid IGP^ This ; m 

arid ovSi? tlae region rirbtind the swamps south west of this town , 
black and vriiite pictures were taken by an RC-10 camera and also 
35 mm color Ektachrome pictures/ This information and that al- 
ready eicistent v?aS used to help both the photointerpretation work 
at ONbrN and the digital analysis at IGP. 

Fig. 14 shoi'Js the result of an early version of INTERPRET 
made during the last part of 1974. Channels 5 and 6 , channel 7 

vjas missing- due to unJ<nown causes > on the scene of Fig. 3 were - 
used. Fig; 14 is a computer output made with multiple print to 
obtain: a gray level presentation. It was thought that this .type 
of presentation would he advantageous for the delineation of the 
different classes. ; . , *■■■’■ 

Fig. 15 shows a color version of this analysis, it was ob- 
tained from computer printouts by coloring each character with 
different colors. Tills figure shows a branch of the Ucayali Ri- 
ver west of Contamana and the swamp area south of it, 

Eiue , violet and black indicates different tones of water 
as seen from Lands at 1. The blade represents highly sedimented - 
water and the other colors, water with different concentrations of 
matter in suspension. 
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FIG. 14 .- Map from an early version of INTERPRET 
Multiple print was used to obtain a 
gray level presentation for the dif- 
ferent classes. 
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Th 0 green color represents a combination of water and vege. 
tation, this was concluded after examining the ground truth and 
considering that the MSS detectors picks up an average of the re- 
flected radiation over an area of aproximately 5,000 tn^ an<5 ther^ 
fore a mixture of water and vegetation is observed at the border 
of lakes and rivers, it is interesting to observe that very na — 
rrow rivers lil<e the one that goes from the swamps toward the 
north also appears green due to this effect on figure 15. 

Orange represents zones where herbaceous vegetation is grow 
ing usually over former meanders or lakes. The dark yellow color 
represents swamp or muddy area. 

Fig. 16 show a result from the present version of INTERPRET 
as it was displayed on the oscilloscope and Fig. 17 a display on 
the teleprinter. In the later picture, special characters design 
ed to produce a gray level presentation were used. 

Fig. 18 shows a color version of the same result, here the 
jungle was also included, an attempt was made to match the classi- 
fication obtained with INTERPRET to that obtained by Malleaux*from 
analysis of aercphotography. Much of the detail obtained with the 
computer was averaged out so that to have a more general presenta- 
tion suitable for this comparison. 


* Forestry Map of Peru 1976. 
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Color map from an early version of INTERPRET, 
The following color code was used 


Black 

Blue and Violet 

Orange. 

Green 


higly sedimented water 

water with intermediate degree 
of sedimentation 

herbaeeous- vegatation, 

mixture of water and vegetation 
and ■ wet areas 


Dark yellow 


swamps 




fir ‘.W' 'V, 


AREA CON VEGETACION ARBOREA 
POBRE 


C0MP05ICI0N OE LAS 3 AREAS 
ANTERIORES 


FIG. 16 - OSCILLOSCOPE MAPS FROM PRESENT 
VERSION OF INTERPRET 















Fig. 18 Color* map from present version of INTERPRET. The 

follov7ing color code was used: 


Yellow 

Alluvial forest type III 

Green 

Mixture of vegetation 
wet areas or alluvial 
II 

and water, 
foi’est type 

Orangs 

Herbaceous vegetation 
tural land 

or agrlcul- 

Blue 

Water 


.Pinlc 

Swamps 


Wiite 

Unrecognized pixels 
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Pig. 19 shows an extract from the "Mapa Porestal del Peru” 
published by UNA at a 1:1,000,000 scale. A comparison with fi- 
gure 18 shows that there are some obvious differences on the 

drainage patterns due to the fact that the aerophotography used - 
as a base for the forestry map is over 20 years old while the 

Landsat imagery is recent, October 1972. It is interesting to 

observe that the lake east to the swamp area, zone I, has almost 
dryed out and an oxbow lake has been formed south Contamana given 
a clear indication of the changes that are continously occuring in 
the meander valley of the Ucayali river. 

It can be seen that in spite of the averaging performed on 
the computer results to obtain figure 18, there is still a great 
deal of detail in this figure when compared with the forestry map. 

'J*he matching of figure 18 and 19 is satisfactory, the 

yellow color represents what has been defined as an alluvial for 
est type III on the forestry map. 

The wet areas or mixture of vegetation and water are repre 
sciited by green color in figure 18. It appears that the large 
arta colored green at the left hand bottom corner of this figure 
mu/ be also a representation of an Alluvial forest type II,, this 
shall be investigated during the continuation of this reseach. 

The orange color represents areas with herbaceous vegeta — 
tion or also areas where the action of man' has talcen place. Blue 
color represents water and pink represents the swamp area-.'/. White 
are those samples that were not recognized by the computer. 



<y/BA-m 
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Ag AGRICULTURA 
Ap AREA PANTANOSA 


BArl BO^UE ALUVIAL I 

I 

BA-11 BOSQUE ALUVIAL II 


BArlll BOSQUE ALUVIAL III 


BO I BOSQUE CO LIN A I 


BC-II BOSQUE CO LIN A II 


BC-III BOSQUE COLINA HI 


BPI BOSQUE PROTECCION I 


EXTRACT FROM "MAPA FORES7AL DEL PERU 1976 
SCALE 1:1,000,000 

NATIONAL AGRARIAN UNIVERSITY 





When the areas south of Contamana are observed from an air 
plane they show an Incipient development of agriculture which co£ 
relates very well with the orange area to the right of figure 18 
South of the. island and, extending toward the South. This is part 
of the area that is classified as agricultural land in the fores- 
try map. 

INTERPRET is also being tested on the study of the city of 
Tarapoto trying to determine the maximum amount of information - 
that can be obtained from a Landsat imagery. Fig. 2:0 shows a pre- 
liminary result obtained on the oscilloscope from zone II. In 

Fig. l3 the city has been clearly delineated along with tne main 
roads that converge on it. 

At present there is not a recent map of the city and its 
• surroundings,, and it is probable that the present development of 
the agricultural areas close to Tarapoto have not been assesed - 
yet. Classification maps from Landsat could be very valuable for 
the determination of the present state of development of cities - 
like Tarapoto. This work with INTERPRET will continue in order 
to determine its actual potential as a tool for demographic stu- 
dies. 

VIII. Conclusions 

The use of clustering methods permit the development of 
relatively fast classification algorithms that could be implemenjt 
ed in an inexpensive computer system with limited amount of iremory. 




FI(j’. 2U - Map of Tarapoto. Brighter points 
I’epi’e.sent urban areas. Medium le, 
vei poircs represent unpaved roods 
Weak level points represent uater. 
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Tlie analysis of CCTs by means of -these teehnigues can provide a 
great deal of detail permitting ‘the use of the maximum resolution 
of Lands at imagery. 

Under the supervision of an specialized user the programs 
developed can produce satisfactory thematic maps in a relatively 
short time if an appropriate recording system is used. 

In the course of this work it has been detected potential 
cases in which 'the use of other techniques for classification us- 
ing a gaussian approximation for the distribution functions can 
'be used with advantage. In general a system that uses several - 
different methods of classification is desirable. 

For the study of jungle areas*, channels 5 and 7 (or 6) can 
provide enough information to delineate drainage patterns , swamp 
and wet areas and make a reasonable broad classification of for- 
est types. 

The results obtained by means cf INTERPRET were found sati^ 
factory, by direct observation of the area from an airplane, aero, 
photography, and when compare with the Forestry Map published by 
UNA. 

Computerized methods sustained by special purpose hardware 
can accelerate the process of evaluation of natural resources and 
represent the most logical way to deal with multiespeotral im- 
agery. 


original page is 
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IX. Future Work 

At present techniques based on the maximum likelyhood ra- 
tio are being developed and will be added to INTERPRET, 

Under the axispices of the Canadian International Develop- 

* 

ment Agency the present computer system will supplemented as 
shown in Fig. 21- This system will permit to produce thematic 
maps on film so that to permit enlargements to an scale appropri- 
ate for practical purposes. 

Techniques for geometric correction, scale and earth rota- 
tion skew, shall be develped when the envisioned additions to the 
present computer system materialize. 

Techniques based on linear and nonlinear adaptive algo 

rithms are also being study as additions to INTERPRET. 


OF POOR QUALITY 
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ABSTRACT 


Landsat multiespectral images and aerophotography was used to 
clasaiify the different jungle types in the region close to the 
cities of Contaraana and Tarap^to. The attached map shows the 
results of this work. 


